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Abstract 


Studies of the spin, parity and tensor couplings of the Higgs boson in the H — ZZ* — 4€, 
H  WW* — evuv and H — yy decay processes at the LHC are presented. The invest- 
igations are based on 25 fb ! of pp collision data collected by the ATLAS experiment at 
Vs = 7 TeV and ys = 8 TeV. The Standard Model (SM) Higgs boson hypothesis, corres- 
ponding to the quantum numbers J^ = 0*, is tested against several alternative spin scenarios, 
including non-SM spin-0 and spin-2 models with universal and non-universal couplings to 
fermions and vector bosons. All tested alternative models are excluded in favour of the SM 
Higgs boson hypothesis at more than 99.9% confidence level. Using the H — ZZ* — 4t and 
H  WW* — evyv decays, the tensor structure of the interaction between the spin-0 boson 
and the SM vector bosons is also investigated. The observed distributions of variables sens- 
itive to the non-SM tensor couplings are compatible with the SM predictions and constraints 
on the non-SM couplings are derived. 
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1 Introduction 


The discovery of a Higgs boson by the ATLAS [1] and CMS [2] experiments at the Large Hadron Collider 
(LHC) at CERN marked the beginning of a new era of experimental studies of the properties of this 
new particle. In the Standard Model (SM), the Higgs boson is a CP-even scalar particle, J“? = 0**. 1 
Theories of physics beyond the SM (BSM) often require an extended Higgs sector featuring several 
neutral Higgs bosons. Such cases may include CP-mixing in the Higgs boson interactions, which could 
result in observable differences in the kinematics of final-state particles produced in their decays. A 
review of the phenomenology in the determination of Higgs boson spin and CP properties can be found 
in Ref. [3] and references therein. 


Previous determinations of the Higgs boson spin and CP quantum numbers by the ATLAS and CMS 
Collaborations are reported in Refs. [4] and [5]. Results on the same subject have also been published by 
the DO and CDF Collaborations in Ref. [6]. All these studies indicate the compatibility of the spin and 
CP properties of the observed Higgs boson with the SM predictions. The ATLAS measurement excluded 
several alternative spin and parity hypotheses in favour of the quantum numbers predicted by the SM. In 
addition to the exclusion of several non-SM spin hypotheses, the CMS measurement probed the tensor 
structure of the Higgs boson decay to SM vector bosons in the spin-0 scenario. This paper complements 
the previous ATLAS study of the Higgs boson spin and parity. The new study takes advantage of im- 
provements to the analysis strategy and to the modelling used to describe alternative spin hypotheses, and 
includes studies on CP-mixing for the spin-0 scenario. The improved theoretical framework is based on 
the Higgs boson characterisation model described in Refs. [3, 7]. 


The study of the spin and parity properties of the Higgs boson presented in this paper is based on the 
H — yy, H — ZZ — 4t and H —^ WW* — evyv decay channels and their combination. The 
H — WW* > evpy analysis is described in detail in a separate publication [8]. These analyses are based 
on 4.5 fb-1 and 20.3 fb^! of pp collision data collected by the ATLAS experiment at centre-of-mass 
energies of 7 TeV and 8 TeV, respectively. For the H —^ WW* — evyy studies only the data collected 
at a centre-of-mass energy of 8 TeV are used. The SM hypothesis J^ = 0* is compared to alternative 
spin-0 models: a pseudoscalar boson J^ = 0^ and a BSM scalar boson J? = 07 [9, 10], which describes 
the interaction of the Higgs boson with the SM vector bosons with higher-dimension operators discussed 
in Section 3.1. Graviton-like tensor models with J? = 2* with universal and non-universal couplings [3, 
7] are also considered. In these tests of fixed spin and parity hypotheses it is assumed that the resonance 
decay involves only one CP eigenstate. 


In addition to the fixed spin and parity hypothesis tests, the possible presence of BSM terms in the Lag- 
rangian describing the HVV vertex? of the spin-0 resonance is also investigated. The HVV interaction 
is described in terms of an effective Lagrangian that contains the SM interaction and BSM CP-odd and 
CP-even terms [3, 7]. The relative fractions of the CP-odd and CP-even BSM contributions to the ob- 
served Higgs boson decays are constrained, and limits on the corresponding BSM tensor couplings are 
derived. 


This paper is organised as follows. In Section 2 the ATLAS detector is described. In Section 3 the 
theoretical framework used to derive the spin and parity models, as well as the parameterisation used 
to describe the HVV coupling tensor structure, are discussed. In Section 4, the choice of Monte Carlo 
generators for the simulation of signal and backgrounds is described. The analyses of fixed spin and parity 


! [n the following, for brevity, only the J? label is used to indicate the spin and CP quantum numbers. 
? [n this paper the symbol V is used to describe a massive SM vector boson, namely either a W or a Z boson. 


hypotheses for the three decay channels and their combination are presented in Section 5. Individual and 
combined studies of the tensor structure of the HVV interaction are presented in Section 6. Concluding 
remarks are given in Section 7. 


2 The ATLAS detector 


The ATLAS detector is described in detail in Ref. [11]. ATLAS is a multi-purpose detector with a forward- 
backward symmetric cylindrical geometry. It uses a right-handed coordinate system with its origin at the 
nominal interaction point (IP) in the centre of the detector and the z-axis along the beam pipe. The 
x-axis points from the IP to the centre of the LHC ring, and the y-axis points upward. Cylindrical co- 
ordinates (r, ó) are used in the transverse plane, ó being the azimuthal angle around the beam pipe. The 
pseudorapidity is defined as 7 = — Intan(0/2), where 0 is the polar angle. 


At small radii from the beamline, the inner detector (ID), immersed in a 2 T magnetic field produced by 
a thin superconducting solenoid located in front of the calorimeter, is made up of fine-granularity pixel 
and microstrip detectors. These silicon-based detectors cover the range |n| < 2.5. A gas-filled straw-tube 
transition-radiation tracker (TRT) complements the silicon tracker at larger radii and also provides elec- 
tron identification based on transition radiation. The electromagnetic (EM) calorimeter is a lead/liquid- 
argon sampling calorimeter with an accordion geometry. The EM calorimeter is divided into a barrel 
section covering |7| < 1.475 and two end-cap sections covering 1.375 < |n| < 3.2. For |n| < 2.5 it is 
divided into three layers in depth, which are finely segmented in 7 and ¢. An additional thin presampler 
layer, covering |7| « 1.8, is used to correct for fluctuations in energy losses of particles before they reach 
the calorimeter. Hadronic calorimetry in the region [| < 1.7 uses steel absorbers and scintillator tiles 
as the active medium. Liquid argon with copper absorbers is used in the hadronic end-cap calorimeters, 
which cover the region 1.5 < |ņ| < 3.2. A forward calorimeter using copper or tungsten absorbers with 
liquid argon completes the calorimeter coverage up to |7| = 4.9. The muon spectrometer (MS) measures 
the deflection of muon trajectories with |7| < 2.7, using three stations of precision drift tubes, with cath- 
ode strip chambers in the innermost layer for [i| > 2.0. The deflection is provided by a toroidal magnetic 
field with an integral of approximately 3 Tm and 6 Tm in the central and end-cap regions of the ATLAS 
detector, respectively. The muon spectrometer is also instrumented with dedicated trigger chambers, the 
resistive-plate chambers in the barrel and thin-gap chambers in the end-cap, covering |7| < 2.4. 


3 Theoretical models 


In this section, the theoretical framework for the measurements of the spin and parity of the resonance 
is discussed. An effective field theory (EFT) approach is adopted to describe the interaction between 
the resonance and the SM vector bosons, following the Higgs boson characterisation model described in 
Refs. [3, 7]. Three possible BSM scenarios for the spin and parity of the boson are considered: 


e the observed resonance is a spin-2 particle, 
e the observed resonance is a pure BSM spin-0 CP-even or CP-odd Higgs boson, 


e the observed resonance is a mixture of the SM spin-0 state and a BSM spin-0 CP-even or CP-odd 
state. 


The third case would imply CP-violation in the Higgs sector. In the case of CP mixing, the Higgs boson 
would be a mass eigenstate, but not a CP eigenstate. In all cases, only one resonance with a mass of 
about 125 GeV is considered. It is also assumed that the total width of the resonance is small compared 
to the typical experimental resolution of the ATLAS detector (of the order of 1—2 GeV in the four-lepton 
and yy final states, as documented in Ref. [12]). Interference effects between the BSM signals and SM 
backgrounds are neglected. 


The EFT approach, used by the Higgs boson characterisation model, is only valid up to a certain energy 
scale, A. The models described in Ref. [7] assume that the resonance structure corresponds to one new 
boson (X(J^) with J? = 0* or 2*), assuming that any other BSM particle only exists at an energy scale 
larger than ^. The A scale is set to 1 TeV to account for the experimental results obtained at the LHC and 
previous collider experiments, which do not show any evidence of new physics at lower energy scales. 


The case where the observed resonance has J^ = 1* is not studied in this paper. The H — yy decay is 
forbidden by the Landau-Yang theorem [13, 14] for a spin-1 particle. Moreover, the spin-1 hypothesis 
was already studied in the previous ATLAS publication [4] in the H —^ ZZ* — 4t and H — WW* — evyuv 
decays and excluded at a more than 99% confidence level. 


3.1 The spin-0 hypothesis 


In the spin-0 hypothesis, models with fixed spin and parity, and models with mixed SM spin-0 and BSM 
spin-0 CP-even and CP-odd contributions are considered. In Ref. [7], the spin-O particle interaction with 
pairs of W or Z bosons is given through the following interaction Lagrangian: 


LY = {cos(a)ksm | jauzzZ,Z" + guwwW; W `“ | 


_ 11 cos(a)KyzzZyyZ” + sin(a)kazzZ,y Z^ | (1) 


-44 [eosto»uww Wa, W^" + sin(a)kaww WW" |] Xo. 


Here V" represents the vector-boson field (V = Z,W*), the V" are the reduced field tensors and the 
dual tensor is defined as V!" = Seto Vee. The symbol A denotes the EFT energy scale. The symbols 
KsM. Kyy and Kayy denote the coupling constants corresponding to the interaction of the SM, BSM 
CP-even or BSM CP-odd spin-0 particle, represented by the Xo field, with ZZ or WW pairs. To ensure 
that the Lagrangian terms are Hermitian, these couplings are assumed to be real. The mixing angle a 
allows for production of CP-mixed states and implies CP-violation for a + 0 and a + a, provided the 
corresponding coupling constants are non-vanishing. The SM couplings, gyyy, are proportional to the 
square of the vector boson masses: gyyy « mÇ. Other higher-order operators described in Ref. [7], 
namely the derivative operators, are not included in Eq. (1) and have been neglected in this analysis since 
they induce modifications of the discriminant variables well below the sensitivity achievable with the 
available data sample. 


As already mentioned, for the spin-0 studies the SM Higgs boson hypothesis is compared to two altern- 
atives: the CP-odd JP = 0^ and the BSM CP-even J^ = 0; hypotheses. All three models are obtained 
by selecting the corresponding parts of the Lagrangian described in Eq. (1) while setting all other con- 
tributions to zero. The values of the couplings corresponding to the different spin-0 models are listed in 
Table 1. 


The investigation of the tensor structure of the HVV interaction is based on the assumption that the ob- 
served particle has spin zero. Following the parameterisation defined in Eq. (1), scenarios are considered 


J! Model Values of tensor couplings 


KsM_ KHVV KAVV Q 

0* SM Higgs boson 1 0 0 0 

07 BSM spin-0 CP-even 0 1 0 0 
0 BSM spin-0 CP-odd 0 0 1 nm/2 


Table 1: Parameters of the benchmark scenarios for spin-0 boson tensor couplings used in tests (see Eq. (1)) of the 
fixed spin and parity models. 


where only one CP-odd or one CP-even BSM contribution at a time is present in addition to the SM 
contribution. To quantify the presence of BSM contributions in H — ZZ* and H —^ WW* decays, the 
ratios of couplings (Kavy/Ksm) : tan à and Kgyy/Kksw are measured. Here KAyy and Kgyy are defined as 


follows: i Í 
KAyy = zav and Kyyy = T Kuy, (2) 


where v is the vacuum expectation value [15] of the SM Higgs field. 


The mixing parameters (Kayy/Ksm) : tano and Kyyy/kKsm correspond to the ratios of tensor couplings 
ga/g, and g2/gı proposed in the anomalous coupling approach described in Refs. [9, 10]. To compare 
the results obtained in this analysis to other existing studies, the final results are also expressed in terms 
of the effective cross-section fractions (fg2, 62) and (fg4, 6,4) proposed in Refs. [3] and [9, 10]. Further 
details of these conversions are given in Appendix A. 


The BSM terms described in Eq. (1) are also expected to change the relative contributions of the vector- 
boson fusion (VBF) and vector-boson associated production (VH) processes with respect to the gluon- 
fusion (ggF) production process, which is predicted to be the main production mode for the SM Higgs 
boson at the LHC. For large values of the BSM couplings, at the LHC energies, the VBF production 
mode can have a cross section that is comparable to the ggF process [16]. This study uses only kinematic 
properties of particles from H — VV* decays to derive information on the CP nature of the Higgs boson. 
The use of the signal rate information for different production modes, in the context of the EFT analysis, 
may increase the sensitivity to the BSM couplings at the cost of a loss in generality. For example the ratio 
of the VBF and VH production modes with respect to the ggF one can be changed by a large amount for 
non-vanishing values of the BSM couplings. In the studies presented in this paper the predictions of the 
signal rates are not used to constrain the BSM couplings. 


As described in Section 6.2, only events with no reconstructed jets (the 0-jet category) are used in the H > 
WW* — evpy analysis for the studies of the tensor structure; hence this analysis has little sensitivity to 
the VBF production mode. The H — ZZ* — 4t analysis also has little sensitivity to this production mode 
since it is mainly based on variables related to the four-lepton kinematics. The Boosted Decision Tree 
(BDT) algorithm [17] used to discriminate signals from the ZZ* background, described in Sections 5.4 
and 6.3, includes the transverse momentum of the four-lepton system and is trained on simulated samples 
of ggF-produced signals. An enhancement of the VBF production mode would improve the separation 
between background and signal since it predicts larger values of the transverse momentum spectrum for 
events produced via VBF than via ggF [3]. 


3.2 The spin-2 hypothesis 


In the Higgs boson characterisation model [7], the description of the interaction of a spin-2 particle with 
fermions and vector bosons is described by the following Lagrangian: 


l 
£= |2 TAXI TAXI. (3) 
v f 


The spin-2 tensor field X^" is chosen to interact with the energy-momentum tensors, TY and T of 
any vector boson V and fermion f, as inspired by gravitation theories. The strength of each interaction 
is determined by the couplings xy and «y. In the simplest formulation, all couplings are equal. This 
scenario is referred to as universal couplings (UC), while scenarios with different values of the couplings 
are referred to as non-universal couplings (non-UC). In the UC scenario, the production of a spin-2 
particle in pp collisions is expected to be dominated by QCD processes, with negligible contributions 
from electroweak (EW) processes (i.e. from processes involving EW boson propagators). Simulation 
studies based on MapGRAPH5 aMCGNLO [16] , which implements the Lagrangian described in Eq. (3), 
predict for the production cross section in the UC scenario oEw/7aQcp = 3 x 1074. These studies also 
show that EW production of the spin-2 resonance would occur mainly in association with a massive 
EW boson (WX, ZX). Present observations do not show a dominant VH production mechanism, hence 
suggesting that agw is significantly smaller than o-ocp. This paper considers only QCD production for 
all the spin-2 benchmark scenarios. 


The UC models predict a branching ratio of about 5% to photon pairs and negligible branching ratios 
to massive EW gauge boson pairs, WW“ and ZZ". This prediction is disfavoured by the experimental 
measurements [18—20] and therefore the equality between all couplings « cannot hold. In the benchmark 
scenarios studied in this paper, each of the couplings kw, xz, and xy is assumed to be independent of 
all the other couplings. In the following, the UC scenario only refers to kg = kç, without implying the 
equality for the other « values. 


The simplest QCD production processes, gg — X and qq — X (where q refers to light quarks), yield 
different polarisations for the spin-2 particle X, and hence different angular distributions of its decay 
products. These mechanisms are considered in the model of a graviton-like tensor with minimal coup- 
lings proposed in Refs. [9, 10], which has been studied experimentally in Ref. [4]. The EFT Lagrangian, 
however, also allows for more complex processes with emission of one or more additional partons. For 
instance, processes with one-parton emission, like gg — qX and qg — 4X, can produce a spin-2 state 
through either a qqX or a ggX vertex. When two partons are emitted, as in gg — qqX or qq — qqX, 
the spin-2 production may occur through qqX or ggX vertices, respectively, such that the polarisation of 
X is not uniquely determined by the initial state. Moreover, the EFT also allows for four-leg vertices 
like qqgX. These additional diagrams effectively change the polarisation of the particle X, compared to 
what is assumed by the model in Refs. [9, 10]. As a consequence, the angular distributions of the decay 
products become harder to separate from those expected for a scalar resonance. 


The QCD production of a spin-2 particle is driven by the values of the couplings Kg, Kq. Presently, there 
are no experimental constraints on the ratio x,/k, from observed decay modes, since the separation of jets 
initiated by gluons or by light quarks is experimentally difficult and has not yet been attempted in Higgs 
boson studies. The ratio «;/«; can thus be regarded as a free parameter. When <ç 7 Kg, the spin-2 model 
predicts an enhancement of the tail of the distribution of the transverse momentum, p , of the spin-2 
particle. Such a high-pi tail is not present for the x, = x, (UC) case. As stated before, however, the EFTs 


are valid only up to some energy scale, ^. At higher energies, new physics phenomena are expected to 
enter to regularise the anomalous ultra-violet behaviour. 


In the present analysis, a selection jm « 300 GeV is applied when investigating non-UC scenarios, 
Kg £ Kg. In addition, for the non-UC scenarios, analyses using a tighter selection pr « 125 GeV are also 
performed. This is a conservative choice for the pt selection, as the EFT must describe the physics at 
least up to the mass of the observed resonance. It has been verified that the choice of the Pr selection does 
not affect the results for the UC scenario. Even assuming the n « 300 GeV selection, some choices of 
Kj [Kg produce high- Pr tails incompatible with the observed differential distribution reported in Refs. [21, 
22]. For this reason the investigated range of the kg/kg ratio is limited to between zero and two. The spin-2 
scenarios considered in this study are presented in Table 2. The x, = Kọ model is referred to hereafter as 
the UC scenario. The x, = 0 case implies a negligible coupling to light quarks, whereas the x, = 2x; case 
is an alternative scenario with an enhanced coupling to quarks. 


Values of spin-2 quark and gluon couplings px selections (GeV) 


Kg = Kg Universal couplings _ _ 
Kg = 0 Low light-quark fraction < 300 < 125 
Kq = 2Kg Low gluon fraction < 300 < 125 


Table 2: Choices of the couplings to quarks x, and to gluons x, studied for the spin-2 benchmark scenarios. The 
values of the selection criteria applied to the transverse momentum Pr of the spin-2 resonance are also shown. For 
the UC scenario no pr selection is applied. 


4 Data and simulated samples 


The data presented in this paper were recorded by the ATLAS detector during the 2012 LHC run with 
proton-proton collisions at a centre-of-mass energy of 8 TeV, and correspond to an integrated luminosity 
of 20.3 fb^!. For the H > yy and H — ZZ* — 4 channels, the data collected in 2011 at a centre- 
of-mass energy of 7 TeV corresponding to an integrated luminosity of 4.5 fb^!, are also used. Data 
quality requirements are applied to reject events recorded when the relevant detector components were 
not operating correctly. More than 90% of the recorded luminosity is used in these studies. The trigger 
requirements used to collect the data analysed in this paper are the same as those described in previous 
publications [18-20]. They are only briefly recalled in the following sections. 


The Monte Carlo (MC) samples for the backgrounds and for the SM Higgs boson signal are the same as 
those used for the analyses described in Refs. [18-20], whereas new non-SM signal samples have been 
simulated. An overview of the signal samples is given in Section 4.1. 


The effects of the underlying event and of additional minimum-bias interactions occurring in the same or 
neighbouring bunch crossings, referred to as pile-up in the following, are modelled with Pytuia 8 [23]. 
The ATLAS detector response is simulated [24] using either Grant 4 [25] alone or combined with a 
parameterised Grant 4-based calorimeter simulation [26]. 


4.1 SM Higgs boson and BSM signal samples 


The SM Higgs boson ggF production for all analyses is modelled using the PowHEc-Box [27] generator 
at next-to-leading order (NLO), interfaced to Pvruri4 8 for parton showering and hadronisation and to 
simulate multi-parton interactions. To improve the modelling of the SM Higgs boson pr, a reweighting 
procedure is applied. This procedure applies a weight depending on the pr of the Higgs boson to each 
event. The weights are chosen in order to reproduce the prediction of the next-to-next-to-leading-order 
(NNLO) and next-to-next-to-leading-logarithms (NNLL) dynamic-scale calculation given by the HREs2.1 
program [28, 29]. 


For the H — yy analysis, the signal samples are generated at several values of the Higgs boson mass 
mg around 125 GeV. The samples are used to obtain a parameterisation of the signal yields and of the 
invariant mass distribution of the two-photon system as continuous functions of my (both inclusively and 
for each category in the analysis, as described in Section 5.2). The spin-2 samples are generated using 
the MaApGRAPH5, aMCGNLO [16] program with LO accuracy for zero, one, and two additional partons, 
and with subsequent matching of the matrix-element calculation with a model of the parton shower, 
underlying event and hadronisation, using PvrHiA 6 [30]. 


In the H — ZZ* — 4€ analysis the signal samples representing the production and decay of Higgs 
bosons with spin-0 and different parities are generated as follows. The SM Higgs boson production via 
gluon fusion at the mass my = 125.5 GeV is simulated using the Pownsc-Box generator. For the non- 
SM signals, the decays of the generated Higgs bosons are simulated, according to the Higgs boson parity 
assumptions, using the JHU [9, 10] MC generator at leading order (LO). The spin-2 samples are generated 
using the MapGrapHS_aMC@NLO MC generator, as for the H — yy analysis. 


For the H —^ WW* — evpy analysis, the SM Higgs boson signal is generated at my = 125 GeV using 
the PownEc-Box Monte Carlo generator. The spin-O BSM signal samples are generated using Map- 
GRAPHS aMCGNLO. The signal samples representing the production and decay of Higgs bosons with 
spin-2 are generated using the MapGrapHS_aMC@NLO MC generator, as for the H — yy analysis. 


For studies of the tensor structure of the HVV decay, all simulated signal samples are obtained by using 
the matrix element (ME) reweighting method applied, as explained in the following, to a sample generated 
with non-zero values of the BSM couplings. The reweighting procedure is validated against samples 
produced at different values of the couplings, to ensure that the distributions of the CP-sensitive final- 
state observables and of their correlations are reproduced correctly. For the H — ZZ* — 4€ analysis, 
the MC production is only performed for one set of tensor couplings: gı = 1, go =1+ig4 = 1+ i. 
All other configurations of couplings are obtained by reweighting this sample at generator level. The 
ratios of the corresponding squares of ME values calculated at LO are used as weights. To calculate these 
ME values, the JHUGenME [10] program is used. In the H — WW*  evgv analysis, only one MC 
sample is generated, using MapGrapHS_aMC@NLO with parameters ksm = 1, Kaww = 2, Kuww = 2, 
cos(a) = 0.3, and all other samples are obtained from it by reweighting the events on the basis of the ME 
amplitudes. 


In all the analyses presented in this paper, the mass of the Higgs boson is fixed to 125.4 GeV [12]. 


4.2 Background samples 


The MC simulated samples for the backgrounds, as well as for the determinations of the corresponding 
cross sections, are the same as those adopted in Refs. [18-20]. In the H — yy analysis, the background 
is dominated by prompt yy events, with smaller contributions from y—jet events. For the H — ZZ* — 4€ 
analysis, the major background is the non-resonant ZZ* process, with minor contributions from the tf and 
Z+jets processes. For the H —^ WW* — evygv analysis, the dominant backgrounds are non-resonant W 
boson pair (WW) production, tf and single-top-quark production, and the Z/y* process followed by the 
decay to Tr final states. 


5 Tests of fixed spin and parity hypotheses 


The H — yy and H — ZZ* — 4t analyses are improved with respect to the previous ATLAS publication 
of Ref. [4]. These analyses are described in some detail in the following subsections. The spin and 
parity analysis in the H — WW* — evyv channel has also been improved, as discussed in detail in a 
separate publication [8]. In the following, only a brief overview of this analysis is given. The expected 
and observed results of the individual channels and of their combination are presented in Section 5.5. 


5.1 Statistical treatment 


The analyses rely on discriminant observables chosen to be sensitive to the spin and parity of the signal. 


A likelihood function, (data | J^, u, 6), that depends on the spin-parity assumption of the signal is con- 
structed as a product of conditional probabilities over binned distributions of the discriminant observables 
in each channel: 


Nehann. Nbins 
2 PY > 2 2 
(data | JP, u, Ê = | | | | P(N; | ijr. (8) + B, (0) - AKO), (4) 
j i 


where u; represents the parameter associated with the signal rate normalised to the SM prediction in each 
channel j.° The symbol 6 represents all nuisance parameters. The likelihood function is a product of 
Poisson distributions P corresponding to the observation of N; ; events in each bin i of the discrimin- 


M 

ant observables, given the expectations for the signal, S 1 6), and for the background, B; FO} Some 

of the nuisance parameters are constrained by auxiliary measurements. Corresponding constraints are 
g > 

represented by the functions A ;(@). 


While the couplings are predicted for the SM Higgs boson, they are not known a priori for the alternative 
hypotheses, defined as P as discussed in Section 3. In order to be insensitive to assumptions on the 
couplings of the non-SM resonance (the alternative hypotheses) to SM particles, the numbers of signal 
events in each channel, for each different LHC centre-of-mass energy and for each tested hypothesis, are 
treated as independent parameters in the likelihood and fitted to the data when deriving results on the spin 


and parity hypotheses. 


3 Here channel can be used to indicate different categories in the same final state when producing results for individual decay 
channels, or different final states when combining them. 


The test statistic g used to distinguish between the two spin-parity hypotheses is based on a ratio of 
profiled likelihoods [31, 32]: 


P & A 
LU sy Bp. Oye) 


— (5) 
Lg ? Hn > a) 


q = log 


where £(J^, ñ JP, à j^) is the maximum-likelihood estimator, evaluated under either the SM Je. = 0* or 


the alternative J A spin-parity hypothesis. The parameters À ye and ó ype represent the values of the signal 
strength and nuisance parameters fitted to the data under each spin and parity hypothesis. The distributions 
of the test statistic for both hypotheses are obtained using ensemble tests of MC pseudo-experiments. For 
each hypothesis test, about 70000 pseudo-experiments were generated. The generation of the pseudo- 
experiments uses the numbers of signal and background events in each channel obtained from maximum- 
likelihood fits to data. In the fits of each pseudo-experiment, these and all other nuisance parameters are 
profiled, i.e. fitted to the value that maximises the likelihood for each value of the parameter of interest. 
When generating the distributions of the test statistic for a given spin-parity hypothesis, the expectation 
values of the signal strengths are fixed to those obtained in the fit to the data under the same spin-parity 
assumption. The distributions of g are used to determine the corresponding p-values [HON = p™ and 
p) = pš. For a tested hypothesis J: the observed (expected) p-values are obtained by integrating 
the corresponding distributions of the test statistic above the observed value of 4 (above the median of the 
JL, Q distribution). When the measured data are in agreement with the tested hypothesis, the observed 


SM 
value of d is distributed such that all p-values are equally probable. 


Very small values of the integral of the distribution of the test statistic for the J a hypothesis, correspond- 
ing to large values of g, are interpreted as the data being in disagreement with the tested hypothesis in 


favour of the SM hypothesis. 


The exclusion of the alternative J A hypothesis in favour of the SM das hypothesis is evaluated in terms 
of the modified confidence level CL,(J7.), defined as [33]: 


alt 


p) 


CL(Jl) = ——L—. 
uii mm 


(6) 


5.2 Spin analysis in the H — yy channel 


The analysis in the H — yy channel is sensitive to a possible spin-2 state. Since the spin-2 models 
investigated in the present paper are different from those assumed in Ref. [4], the analysis has been 
redesigned, to improve its sensitivity to the new models. 


The selection of H — yy candidate events is based on the procedure of other recent ATLAS H yy 
analyses (see for example Ref. [20]). Events are selected if they satisfy a diphoton trigger criterion 
requiring loose photon identification, with transverse momentum pr thresholds of 35 GeV and 25 GeV 
for the photon with the highest (y1) and second-highest (y2) pr, respectively. During the offline selection 
two photons are further required to be in a fiducial pseudorapidity region, defined by |7”| < 2.37, where 
the barrel/end-cap transition region 1.37 < |ņ”| < 1.56 is excluded. The transverse momentum of the 
photons must satisfy p > 0.35 - m, and pr > 0.25 - my, and only events with a diphoton invariant 
mass myy between 105 GeV and 160 GeV are retained. For the events passing this selection, a further 
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requirement is applied on the diphoton transverse momentum, p « 300 GeV, motivated by the assumed 
validity limit of the spin-2 EFT model, as explained in Section 3. After this selection, 17 220 events are 
left at a centre-of-mass energy ys = 7 TeV and 94540 events at s = 8 TeV. 


Kinematic variables sensitive to the spin of the resonance are the diphoton transverse momentum p and 
the production angle of the two photons, measured in the Collins-Soper frame [34]: 


[cos ^| =  |sinh(Ag?)| 2p Pr (7) 
yy 2 my 
1+ (p; [myy) YY 
where Aj?" is the separation in pseudorapidity of the two photons. 
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Figure 1: Expected distributions of kinematic variables sensitive to the spin of the resonance considered in the 
H Ə yy analysis, (a) transverse momentum of the yy system p? and (b) the production angle of the two photons 
in the Collins-Soper frame | cos 6*|, for a SM Higgs boson and for spin-2 particles with three different choices of 
the QCD couplings. 


The predicted distributions of these variables, for events passing the selection, are shown in Figure 1, for 
a SM Higgs boson and for a spin-2 particle with different QCD couplings. For the kq 7 Kg cases, the 
enhanced high-p7" tail offers the best discrimination, whereas for Kç; = Kg the most sensitive variable is 
|cos 6*]. 


To exploit the signal distribution in both pe and | cos 0*|, the selected events are divided into 11 mutually 
exclusive categories: 10 categories (labelled from C1 to C10) collect events with pe < 125 GeV, divided 
into 10 bins of equal size in | cos 6*|, while the 11" category (labelled C11) groups all events with p 2 
125 GeV. As described in Section 3, for the non-UC spin-2 models the analysis is performed with two 
p. selections, namely p « 300 GeV and p? < 125 GeV: the latter case corresponds to not using the 
11* category. 


The number of signal events above the continuum background can be estimated through a fit to the 
observed my, distribution in each category. The m,, distribution is modelled in each category as the sum 
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of one-dimensional probability density functions (pdf) for signal and background distributions: 


"a p (my) + (nl + n) pe (my) 


fl ny) = : (8) 
[c] [c] [c] 
Np thy + Mya 


where J is the spin hypothesis, ni and nl are the background and the signal yield in category c, and 
f, PRG fi my) are the m,, pdfs for the background and the signal, respectively. The signal pdf 


rod any) is modelled as a weighted sum of a Crystal Ball function, describing the core and the lower 
mass tail, and of a Gaussian component that improves the description of the tail for higher mass values. 
For each category, fll my) Is fitted to the simulated my, distribution of the SM Higgs boson and verified 
to be consistent also with the spin-2 models. The background pdf f PIC is empirically modelled as an 
exponential of a first- or second-degree polynomial. The choice of such a Pace en can induce a 
bias ("spurious signal") in the fitted signal yield, which is accounted for by the term nie . The size of the 
expected bias is determined as described in Refs. [20, 22], and ranges between 0.6 and 4 events, depending 
on the category (with the signal ranging from 15 to more than 100 events). In the statistical analysis, au. 
is constrained for each category by multiplying the likelihood function by a Gaussian function centred at 
zero and with a width determined by the size of the expected bias. 


Defining ns as the total signal yield (summed over all categories), the expected fraction of signal events 
[c] 


[c] 


belonging to each category, oi zs depends on the spin hypothesis J. The values of P,- extracted 
n 


from the data can be compared to their expected values for each spin hypothesis, as shown in Figure 2 for 
the data collected at ys = 8 TeV. 
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Figure 2: Observed signal fraction per category for the H — yy analysis, and comparison to expected values for a 
SM Higgs boson and for a spin-2 particle with different choices of QCD couplings. (a) the 11 categories described 
in the text are displayed, corresponding to the p" < 300 GeV selection; (b) the high-p}” category is discarded 
and the signal fractions are renormalised over the 10 remaining categories, corresponding to the p « 125 GeV 
selection. 


For the non-UC scenario the 11^^ (high- pr 7) category provides strong discrimination power against the 
non-SM hypothesis, as visible in Figure 2(a). 


To discriminate between the SM spin-0 Uu = 0*) and alternative spin-2 hypotheses (JP ), two likelihood 


alt’? 
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functions £ Jb £ jh are built, following the general approach described in Eq. (4): 


-n£;-2) (nl; cl + ng qi enl) 


i 
c ec[c] 


-» In [ni fi om) + (ns ql + nicl ) fmt) )| (9) 


bias) 


where >’, runs over all categories and >’ ,<;-] runs over all events in category c. The total signal yield ns is 
a free parameter in the likelihood model. The spin hypothesis being tested enters the likelihood function 
through the fractions of signal per category, b 


Several systematic uncertainties enter this model. They are implemented for each spin hypothesis as 
nuisance parameters, 0j, constrained by multiplicative Gaussian terms in the likelihood function (not 
included in Eq. (9) for simplicity). 


The signal fractions, qi, for the SM Higgs boson are affected by uncertainties on the pr spectrum of 
the resonance and on the size of the interference between the resonance and continuum production. The 
former is computed as described in Ref. [20]. The relative impact on the signal fractions is less than +1% 
for categories 1 to 8 (qp < 125 GeV and | cos @*| < 0.8), and becomes as large as +13% for categories 
10 and 11. The correction for the interference is evaluated according to Refs. [35, 36]. The systematic 
uncertainty is conservatively assumed to equal the correction itself, and its relative impact ranges between 
x 0.146 and «1.8406. 


No systematic uncertainty is assigned to the simulated px distribution of the spin-2 models. The effect of 
the interference between the resonance and continuum production is essentially not known, as it depends 
on the width, Iy, of the resonance, which is unknown. The results presented here only hold under the 
assumption of a narrow width for the resonance, such that interference effects can be neglected. 


Additional systematic uncertainties come from the calibration of the photon energy scale and energy 
resolution and affect the signal parameterisation f; [c]. These uncertainties are evaluated as described in 
Ref. [12]. 


5.3 Spin and parity analysis in the H —^ WW* — evyv channel 


The analysis of the spin and parity in the H —^ WW* — evuv channel is described in detail in a separate 
publication [8]. In the following a brief summary is provided. The selection is restricted to events contain- 
ing two charged leptons of different flavour (one electron and one muon). The evuv channel is the most 
sensitive one [19]. The same-flavour channels (evev and uvuv) are not expected to add much in terms 
of sensitivity due to the presence of large backgrounds that cannot be removed without greatly reducing 
the acceptance of the alternative models considered in this analysis. The leading lepton is required to 
have pr > 22 GeV and to match the object reconstructed by the trigger, while the sub-leading lepton 
needs to have pr > 15 GeV. While the spin-0 analyses select only events with no jets in the final state 
(no observed jets with pr > 25 GeV within [| < 2.5 or with pr > 30 GeV within 2.5 < |n| < 4.5), the 
spin-2 analysis enlarges the acceptance by allowing for zero or one jet (selected according to the above 
mentioned criteria). 


The major sources of background after the dilepton selection are Z/y*+jets (Drell-Yan) events, diboson 
(WW, WZ/y*, ZZ/y*), top-quark (tt and single top) production, and W bosons produced in association 
with hadronic jets (W jets), where a jet is misidentified as a lepton. The contribution from misidentified 
leptons is significantly reduced by the requirement of two high-pr isolated leptons. Drell-Yan events 
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are suppressed through requirements on some of the dilepton variables* (pf > 20 GeV, Adee < 2.8), 
while a cut on mee (mee < 80 GeV) targets the WW background. For alternative spin models with non- 
universal couplings, as discussed in Section 3, an additional upper bound is imposed on the Higgs boson 
pr, reconstructed as the transverse component of the vector sum of the momenta of the two charged 
leptons and the missing transverse momentum. Additionally, for events containing one jet, which include 
substantial top-quark and W+jets backgrounds, b-jet and Z — t*t™ vetoes are applied, together with 
transverse mass requirements: the larger of the transverse masses of the two W bosons (each computed 
using the corresponding lepton and the missing transverse momentum) in the event is required to be 
larger than 50 GeV, while the total transverse mass of the WW system (defined with the two leptons and 
the missing transverse momentum) is required to be below 150 GeV. 


Control regions (CRs) are defined for the WW, top-quark and Drell-Yan backgrounds, which are the 
most important ones after the topological selection described above. The CRs are used to normalise 
the background event yields with a fit to the rates observed in data. The simulation is then used to 
transfer these normalisations to the signal region (SR). The W+jets background is estimated entirely from 
data, while non-WW diboson backgrounds are estimated using MC simulation and cross-checked in a 
validation region. 


After the signal region selection, 4730 and 1569 candidate events are found in data in the 0-jet and 1-jet 
categories, respectively. For the latter category, the number decreases to 1567 and 1511 events when 
applying a selection on the Higgs boson pr of less than 300 GeV and less than 125 GeV, respectively. In 
total 218 (77) events are expected from a SM Higgs boson signal in the 0-jet (1-jet) category, while about 
4390 (1413) events are expected for the total background. 


A BDT algorithm is used in both the fixed spin hypothesis tests and the tensor structure analyses. For 
spin-2 studies, the strategy follows the one adopted in Ref. [4], with the main difference being that the 
1-jet channel has been added. Two BDT discriminants are trained to distinguish between the SM hy- 
pothesis and the background (BDT»), and the alternative spin hypothesis and the background (BDT»). 
Both BDTs employ the same variables, namely mee, pit , Adee and mr, which provide the best discrim- 
ination between signal hypotheses and backgrounds, also in the presence of one jet in the final state. All 
background components are used in the trainings. In total, five BDT, trainings are performed for the 
alternative spin hypotheses (one for the spin-2 UC scenario and two for each of the two spin-2 non-UC 
hypotheses corresponding to the different pe selections), plus one training of BDT 9 for the SM Higgs 
boson hypothesis. 


For the spin-0 fixed hypothesis test and HWW tensor structure studies, the first discriminant, BDTo, is 
the same as the one used for the spin-2 analysis, trained to disentangle the SM hypothesis from the back- 
ground. A second BDT discriminant, BDTep, is obtained by training the SM signal versus the alternative 
signal sample (the pure CP-even or CP-odd BSM hypotheses), and then applied to all CP-mixing frac- 
tions. No background component is involved in this case. The variables used for the BDTe trainings are 
mee, Adee, pr and the missing transverse momentum for the CP-even analysis and mee, Adee, Eeeyy and 
Apr for the CP-odd analysis. The training strategy is different from the one used in the spin-2 analysis 
because, while the spin-2 signal is very similar to the background, the spin-O signals are all similar to 
each other, while being different from the main background components. Therefore, in the latter case, 
training the signal hypotheses against each other improves the sensitivity. The resulting BDT variable is 


4 Throughout this section, the following variables are used: pit and me are the transverse momentum and the invariant mass 
of the two-lepton system, respectively, Aó;; is the azimuthal angular difference between the two leptons, mr is the transverse 
mass of the reconstructed Higgs boson decay system, Apr is the absolute value of the difference between the momenta of the 
two leptons and Egey = pi -0.5 p? + 0.5p7"*, where p7/* is the missing transverse momentum. 
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afterwards used in binned likelihood fits to test the data for compatibility with the presence of a SM or 
BSM Higgs boson. 


Several sources of systematic uncertainty are considered, both from experimental and theoretical sources, 
and are described in detail in Ref. [8]. The correlations induced among the different background sources 
by the presence of other processes in the control regions are fully taken into account in the statistical 
procedure. The most important systematic uncertainties are found to be those related to the modelling of 
the WW background, to the estimate of the W+jets background (originating from the data-driven method 
employed) and, for the spin-2 results in particular, to the Z — tr modelling. 


5.4 Spin and parity analysis in the H — ZZ* — 4€ channel 


The reconstruction of physics objects and event selection used for the H — ZZ* — 4€ analysis is identical 
to the one presented in Ref. [12]. The main improvement with respect to the previous ATLAS publication 
of Ref. [4] is the introduction of a BDT discriminant designed to optimise the separation between the 
signal and the most relevant background process. 


Events containing four reconstructed leptons (electrons or muons) in the final state are selected using 
single-lepton and dilepton triggers. The selected events are classified according to their final state: 
4u, 2e2u, 2u2e and 4e, where for the decay modes 2e2u and 2y2e the first pair is defined to be the 
one with the dilepton mass closest to the Z boson mass. Each muon (electron) must satisfy pr > 6 GeV 
(pr > 7 GeV) and be measured in the pseudorapidity range || < 2.7 (\n| < 2.47). Higgs boson can- 
didates are formed by selecting two same-flavour, opposite-charge lepton pairs in an event. The lepton 
with the highest pr in the quadruplet must have pr > 20 GeV, and the leptons with the second- and 
third-highest pr must have pr > 15 GeV and pr >10 GeV, respectively. The lepton pair with the mass 
closest to the Z boson mass is referred to as the leading lepton pair and its invariant mass as m». The 
requirement 50 GeV < mj» < 106 GeV is applied. The other lepton pair is chosen from the remain- 
ing leptons as the pair closest in mass to the Z boson. Its mass, denoted hereafter by m34, must satisfy 
12 GeV < m34 < 115 GeV. Further requirements are made on the impact parameters of the leptons 
relative to the interaction vertex and their isolation in both the tracker and calorimeter. 


The main background process affecting the selection of H — ZZ* — 4€ events is the non-resonant 
production of ZZ* pairs. This background has the same final state as the signal events and hereafter is 
referred to as the irreducible background. It is estimated from simulation and normalised to the expected 
SM cross section calculated at NLO [37, 38]. The reducible sources of background come from Z+jets and 
tt processes, where additional leptons arise due to misidentified jets or heavy-flavour decays. The rate 
and composition of the reducible backgrounds are evaluated using data-driven techniques, separately for 
the two final states with sub-leading muons ££ + uu and those with sub-leading electrons ££ + ee. 


Only events with an invariant mass of the four-lepton system, denoted by mug, satisfying the signal region 
definition 115 GeV < mag «130 GeV are selected. The expected signal and background yields in the 
signal region and the observed events in data are reported in Table 3. 


The choice of production and decay angles used in this analysis is presented in Figure 3, where the 
following definitions are used: 


e 0, and 0; are defined as the angles between final-state leptons with negative charge and the direction 
of flight of their respective Z bosons, in the four-lepton rest frame; 
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SM Signal ZZ tt,Z * jets Total expected Observed 
Vs =7 TeV 
4u 1.02x0.10 — 0.65x0.03 0.14+0.06 1.812x0.12 3 
2u2e  0.47x0.05  0.29x0.02 0.53+0.12 1.29+0.13 1 
2e2u 0.64+0.06 0.4540.02 0.13+0.05 1.22+0.08 2 
4e 0.45x0.04  0.2640.02 0.59+0.12 1.3040.13 2 
Total 2.5840.25 1.6540.09  1.39x0.26 5.62+0.37 8 
Vs = 8 TeV 
4u 5.8140.58  3.36€0.17 0.97+0.18 10.14+0.63 13 
2u2e = 3.0040.30 1.59+0.10 0.52+0.12 5.11+0.34 8 
2e2u = 3.7240.37 2.33+0.11 0.84+0.14 6.89+0.41 9 
4e 2.01x0.29 1.44+0.09 0.52+0.11 4.87+40.32 7 
Total 15.4 x1.5  8.72x0.47  2.85+0.39 27.0 +1.6 37 


Table 3: Expected signal, background and total yields, including their total uncertainties, and observed events in 
data, in the 115 GeV < mag < 130 GeV signal region. The number of expected signal events is given for a SM Higgs 
boson mass of 125.5 GeV. 


e @ is the angle between the decay planes of two lepton pairs (matched to the two Z boson decays) 
expressed in the four-lepton rest frame; 


e @, is the angle between the decay plane of the leading lepton pair and a plane defined by the Zi 
momentum (the Z boson associated with the leading lepton pair) in the four-lepton rest frame and 
the positive direction of the collision axis; 


e 6* is the production angle of the Zi defined in the four-lepton rest frame. 


The final-state observables sensitive to the spin and parity of a boson decaying to ZZ* — 4€ are the two 
production angles 6* and ®, and the three decay angles ®, 6, and 6. In the case of a spin-0 boson, 
the differential production cross section does not depend on the production variables cos(@*) and ®,. It 
should be noted that, as the Higgs boson mass is below 27nz, the shapes of the mass distributions of the 
intermediate Z bosons, mı2 and m34, are sensitive to the spin and parity of the resonance. In Figure 4 the 
distributions of the final-state observables sensitive to the spin and parity of the decaying resonance are 
presented. The distributions are shown for the SM J?” = 0* and JP = 0^ simulated events, as well as for 
ZZ* production and reducible backgrounds in the signal region 115 GeV « mag « 130 GeV. The events 
observed in data are superimposed on each plot. 


Two approaches were pursued to develop the discriminants used to distinguish between different spin 
and parity hypotheses. The first uses the theoretical differential decay rate for the final-state observables 
sensitive to parity to construct a matrix-element-based likelihood ratio analysis (J^-MELA). The second 
approach is based on a BDT. 


For the /"-MELA approach [3, 9], the probability of observing an event with given kinematics can 
be calculated. This probability is corrected for detector acceptance and analysis selection, which are 
obtained from the simulated signal MC samples. The full pdf also includes a term for incorrect pairing 
of the leptons in the 4u and 4e channels. For a given pair of spin-parity hypotheses under test, the final 
discriminant is defined as the ratio of the pdf for a given hypothesis to the sum of the pdfs for both 
hypotheses. 
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Figure 3: Definitions of the angular observables sensitive to the spin and parity of the resonance in the 
X  ZZ* — 4t decay. 


For the BDT approach, a J^ discriminant is formed for each pair of spin-parity states to be tested, 
by training a BDT on the variables of simulated signal events which fall in the signal mass window 
115 GeV < mag < 130 GeV. For the 0* versus O` test, only the parity-sensitive observables ®, 61, 62, 
mio and m34 are used in the BDT training. For the spin-2 test, the production angles 6* and 6, are also 
included. 


Both analyses are complemented with a BDT discriminant designed to separate the signal from the ZZ* 
background. These discriminants are hereafter referred to as BDTzz. For the J^—MELA analysis, the 
BDTzz discriminant is fully equivalent to the one described in Refs. [12, 18]. For the BDT analysis the 
discriminating variables used for the background BDTzz are the invariant mass, pseudorapidity, and 
transverse momentum of the four-lepton system, and a matrix-element-based kinematic discriminant 
Kp defined in Ref. [16]. The results from both methods are obtained from likelihood fits to the two- 
dimensional distributions of the background BDTs and of the spin- and parity-sensitive discriminants. 
In this way, the small correlation between these variables are taken into account in the analyses. The 
distribution of the background discriminant BDTzz versus the J’-MELA discriminant is presented in 
Figure 5 for the SM J^ = 0* signal, the backgrounds, and the data. The projections of this distribution on 
the JP—MELA and the BDTzz variables, for different signal hypotheses, the backgrounds, and the data, 
are shown in Figure 6. In this paper, only results based on the J^—MELA approach are reported. The 
BDT approach was used as a cross-check and produced compatible results. 


Two general types of systematic effects impact the analyses using fixed spin and parity hypotheses: un- 
certainties on discriminant shapes due to experimental effects, and uncertainties on background normal- 
isations from theory uncertainties and data-driven background estimates. The systematic uncertainties on 
the shape are included in the analysis by creating discriminant shapes corresponding to variations of one 
standard deviation in the associated sources of systematic uncertainty. The systematic uncertainties on 
the normalisation are included as additional nuisance parameters in the likelihood. 


The list of sources of systematic uncertainty common to all ATLAS H — ZZ* — 4€ analyses is presented 
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Figure 4: Distributions of some of the final-state observables sensitive to the spin and parity of the resonance in 
the H — ZZ* — 4€ signal region 115 GeV < ma; < 130 GeV for data (points with errors), backgrounds (filled 
histograms) and predictions for two spin hypotheses (SM solid line and alternatives dashed lines). (a), (b) and (c): 
invariant masses mio , m34 and decay cos 6, respectively; (d), (e) and (f): ®, cos0* and Ó, respectively. 


in Ref. [18]. The relative impact of these sources on the final separation for all tested hypotheses is evalu- 
ated and sources affecting the final separation (given in Section 5.5) by less than +0.5% are neglected. 


The main sources of systematic uncertainties are related to the experimental error on the Higgs boson 
mass, the modelling of the irreducible ZZ* background, the uncertainty on the integrated luminosity and 
the experimental uncertainties on the electron and muon reconstruction. The uncertainty on the Higgs 
boson mass affects the final result since it impacts the shapes of the mio, m34, cos 6; and cos 6» variables. 
For the JP-MELA method, the uncertainty on the estimate of the fraction of 4u and 4e candidates with an 
incorrect pairing of leptons is also considered. This uncertainty is derived by comparing the corresponding 
prediction obtained from the Pownec and JHU MC generators for the SM hypothesis. A variation of 
+10% of the incorrect pairing fraction is applied to all spin and parity hypotheses. 


The influence of the main systematic uncertainties on the separation between the SM J^ = 0* and J? = 07 
hypotheses for the J’-MELA analysis is presented in Table 4. The total relative impact of all system- 
atic uncertainties on the separation between the hypotheses (expressed in terms of numbers of standard 
deviations) is estimated to be about +3%. 
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Figure 5: The distributions of the discriminant BDTzz versus the J?’-MELA discriminant for the SM J? = 0* 
Higgs boson and for the backgrounds in the H — ZZ* — 4€ signal region 115 GeV < mag < 130 GeV. 


Source of the systematic uncertainty Relative impact 


Higgs boson mass experimental uncertainty +2% 

ZZ* pdf +0.8% 
Muon momentum scale +0.7% 
Zbb — (tup normalisation +0.6% 
ZZ* scale +0.6% 
Luminosity +0.6% 
e/y resolution model (sampling term) +0.5% 
e/y resolution model (constant term) +0.5% 
Z (tee normalisation x: 0.596 
Fraction of wrongly paired 4¢ candidates x 0.446 


Table 4: Relative impact of the main systematic uncertainties on the expected separation (expressed in terms of 
numbers of standard deviations) between the SM J? = 0* and J? = 0" hypotheses for the H > ZZ* — 4€ 
J"-MELA analysis. 


5.5 Individual and combined results 


The distributions of discriminant variables in data agree with the SM predictions for all three channels, 
and exclusion ranges for alternative spin hypotheses are derived. Some examples of distributions of the 
test statistic d (defined in Section 5.1) used to derive the results are presented in Figure 7. In this figure, 
the observed value is indicated by the vertical solid line and the expected medians by the dashed lines. 
The shaded areas correspond to the integrals of the expected distributions used to compute the p-values 
for the rejection of each hypothesis. The signal strengths per decay channel and per centre-of-mass energy 
are treated as independent parameters in each fit. Their values are compatible with the SM predictions. 
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Figure 6: Distributions of the J’-MELA and of the BDTzz discriminants in the H > ZZ* — 4€ signal region 
115 GeV < mag < 130 GeV for the data (points with errors), the backgrounds (filled histograms), and for predic- 
tions for several spin and parity hypotheses. The SM hypothesis is shown by the solid line while the alternative 
hypotheses are shown by the dashed lines. The signal distributions are normalised to the signal strength fitted in 
data. (a), (b) and (c): J"-MELA discriminants for 0* SM vs 07, 0* SM vs 0; and 0* SM vs 2*, respectively; (d), 
(e) and (f): BDTzz discriminant for 0* SM vs 07, 0* SM vs 0; and 0* SM vs 2^, respectively. 


The results obtained from the fit to the data, expressed in terms of p-values for different tested hypotheses 
and observed CL, for the alternative hypotheses, are summarised in Tables 5 and 6. As shown in Table 5, 
the sensitivity to reject alternative hypotheses is driven by the H — ZZ* — 4€ and the H —^ WW* — evuv 
channels. The H — yy channel has sizeable sensitivity only to spin-2 models where the Pr < 125 GeV 
selection is not applied. In all cases the data prefer the SM hypothesis to the alternative models, with 
the exception of some of the spin-2 models for the H — yy channel. In this case both hypotheses have 
similar observed p-values, but neither of the two is below 10% . 


As summarised in Table 6, the p-values of the combined results for the three channels show good agree- 
ment between the data and the SM hypothesis for all performed tests. All tested alternative hypotheses 
are rejected at a more than 99.9% confidence level (CL) in favour of the SM hypothesis. 


6 Study of CP-mixing and of the HVV interaction tensor structure 


Following the discussion in Section 3, measurements of the HVV interaction tensor couplings KsM, Kavv; 
Kyvy and of the mixing angle o are performed. The measurements consist of fitting the ratios of couplings 
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Figure 7: Examples of distributions of the test statistic g defined in Section 5.1, for the combination of decay 
channels. (a): 0* versus 07; (b): 0* versus p (c): 0* versus the spin-2 model with universal couplings (kg = Ky); 
(d): 0* versus the spin-2 model with x, = 2x, and the pr selection at 125 GeV. The observed values are indicated 
by the vertical solid line and the expected medians by the dashed lines. The shaded areas correspond to the integrals 
of the expected distributions used to compute the p-values for the rejection of each hypothesis. 
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Tested Hypothesis 


2" (Kg = Kg) 

2* (Kg = 0; pr < 300 GeV) 
2*(k; = 0; pr < 125 GeV) 
2*(k, = 2kg; pr < 300 GeV) 
2*(k, = 2kg; pr < 125 GeV) 


2*(k, = 0; pr < 300 GeV) 
2*(k; = 0; pr < 125 GeV) 
2* (Kg = 2kg; pr < 300 GeV) 
2*(Kg = 2kg; pr < 125 GeV) 


Tested Hypothesis ^ Obs. CL, (96) 


2*(k, = 0; pr < 300 GeV) 
2*(k; = 0; pr < 125 GeV) 
2*(k, = 2kg; pr < 300 GeV) 
2*(Kg = 2kg; pr < 125 GeV) 


Table 5: Expected and observed p-values for different spin-parity hypotheses, for each of the three channels H — 
yy, H  ZZ' — 4, and H —^ WW*  evyv. The observed CL, for the alternative hypotheses are reported 
in the last column. The expected and observed p-values and the observed CL, are defined in Section 5.5 and the 
alternative hypotheses are those described in Section 3. 


(Kavv/Ksm) - tan à and Kgyy/ksM to the discriminant observables for the H —n WW* — evuv and H > 
ZZ* — 4 processes and in their combination. In the fitting procedure only one ratio of couplings 
(Kavv/Ksm) : tan a or Kgyy/KsM is considered at a time, while the other one is assumed to be absent. 


6.1 Statistical treatment 


The measurement of the tensor structure of the HVV interaction is based on a profiled likelihood [31, 32] 
that contains the discriminant observables sensitive to the EFT couplings. The signal rates in the different 
channels and for different centre-of-mass energies are treated as independent parameters. "Therefore, 
the global signal normalisation is not used to constrain the EFT couplings. The ratios of the BSM to 
SM couplings, Kyvy/Ksm and (Kavy/Ksm) : tan a, are each separately fit to the discriminant observables in 
data. The test statistic used to derive the confidence intervals on the parameters of interest is q’ = —2 In(A), 
where 4 is the profiled likelihood [31, 32]. The results presented in the following rely on the asymptotic 
approximation [31, 32] for the test statistic. This approximation was cross-checked with Monte Carlo 
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T alt at SM alt Obs. CL, (%) 


ested Hypothesis Dexp, Pi. Bons Pans 
07 1.8. 107° 13.10 4 088 «31.10? «2.6.10? 
2 (0546) 4.3 - 10? 2.9.107^ 0.61 43.10? 1.1 - 107? 


2*(k4 20; pr < 300 GeV) «3.1.10? <3.1-1075 0.52 <3.1-105 «65.10? 
2*(k, = 0; pr < 125 GeV) 3.4- 10? 39.10434 0.71 43-10? 1.5. 107? 
2*(kq = 2kg; pr < 300 GeV) «3.1.10? <3.1-10° 028 «31-10? «43.10? 
2* (kg = 2kg; pr < 125 GeV) 7.8.10? 1.2.1073 080 7.3.10 3.7. 107? 


Table 6: Expected and observed p-values for different spin-parity hypotheses, for the combination of the three 
channels: H — yy, H — ZZ* — At and H —^ WW* — evyv. The observed CL, for the alternative hypothesis is 
reported in the last column. The expected and observed p-values and the observed CL, are defined in Section 5.5. 
The definitions of alternative hypotheses are given in Section 3. 


ensemble tests that confirm its validity in the range of the parameters for which the 95% CL limits are 
derived. 


6.2 Tensor structure analyses in the H >—- WW* — evpy channel 


The H —^ WW* — evuv analysis used to study the spin-O tensor structure is already described in Sec- 
tion 5.3 and detailed in Ref. [8]. Only the 0-jet category is considered and the BDTp and BDT¢» are used 
as discriminant variables in the likelihood defined to measure the spin-O tensor structure couplings. The 
only difference with respect to the spin hypothesis test is that, in this analysis, the BSM spin-0 couplings 
are treated as continuous variables in the test statistic. 


6.3 Tensor structure analyses in the H — ZZ* — 4( channel 


To allow for a cross-check and validation of the obtained results, two different fitting methods based on 
the analytical calculation of the leading-order matrix element of the H — ZZ* — 44 process are used. 


The method of the matrix-element-observable fit is based on modelling the distributions of the final-state 
observables in each bin of coupling ratios using Monte Carlo simulation. Using the Lagrangian defined 
in Eq. (1), which is linear in the coupling constants «sy, Kgyy and Kayy, the differential cross section 
at each point in the phase space can be expressed as a term corresponding to the SM amplitude, plus 
two additional terms, linear and quadratic in the coupling constants. In this way it is possible to define 
two observables for each coupling, the so-called first- and second-order optimal observables, upon which 
the amplitude depends at each point of the phase space. For each event, they contain the full kinematic 
information about the couplings, which can thus be extracted from a fit to their shapes. More details of 
the method can be found in Refs. [39-42]. 


The observables sensitive to the presence and structure of Ksm, Kgyy and kíyy considered in the current 
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analysis are defined as follows: 


2K [ME(ksm Z0; KHyv,KAvyv=0; a=0)*-ME(kyyy Z0; ksM.KAvy 720; a o 
IME(ksm#0; Kv Kavy 70; a0)? 
- IME(kuvv £0; ksm,kavv 20; a0)? 
OxGamvv) = IME(ksm#0; Krrvv.Kavv -0; a0) ° 
2R[ME(ksm Z0; KHVV.KAVV =! =0; Q= 0)* :ME(Xxavy Z0; KSM, KHvv=0; Q=: 2 
IME(ksm#0; Ky Kavy 70; a=0)|2 
IME(kavy #0; ksm.kHvv 70; a=m/2)2 (10) 
IME(ksm#0; Kyvv.kavv=0; e=0)2 ` 


O (Kuvv) = 


O (Kavv, @) = 
O»(Kavy, @) = 


Here ME(ksm, Kyvv; Kavv, œ) denotes the leading-order matrix element of the H — ZZ* — 4€ process. 
These definitions correspond to the first- and second-order optimal observables for a BSM amplitude with 
a three-component structure. 


The observables O; »(kuyy) and Oj 2(Kayv, œ) are used for the Kyyy/Ksm and (KAyy/KsM)-tan o individual 
fits respectively. In order to suppress the ZZ* background, a kinematic BDT discriminant similar to those 
described in Section 5.4 is used as an additional observable in all fits. The BDT training is performed 
independently for each final state using observables with small sensitivity to parity: 74¢, PT,4¢e, Mae, cos(0*) 
and Ó. This BDT discriminant is denoted hereafter by BDT(ZZ). 


To simplify their use in the analysis, all observables defined in Eq. (10) undergo a pdf transformation 
such that each observable becomes normally distributed in the Standard Model case. These transformed 
observables are referred to hereafter as TO; »(kgvyy) and T O1.2(Kavv, œ) respectively. The distributions of 
transformed observables for the Monte Carlo signal samples generated with (Kqyvy/Ksm = 0, +1; ŘKąvy = 
0) and ((KAvy /Ksw)- tan æ = 0, +5; Kgyy = 0) are shown in Figure 8. The contributions of all backgrounds 
considered in this analysis are also included. By construction the T O> observables are sensitive to the 
modulus of the Kyyy/Ksm and (Kayy/KsM) ` tana ratios: their distributions change with the strength of the 
respective coupling. These observables are insensitive to the relative sign of Kyyy and KAyy with respect 
to ksm. The sign sensitivity comes from the TO, observables, which are based on the interference terms: 
their distributions feature pronounced sign-dependent asymmetries. It was also found that the observables 
T Oi(Kgyy) and TO»(Kgyy) are linearly correlated. To maximise the population of analysis histograms 
with currently available Monte Carlo event samples, it is desirable to reduce this correlation. This is 
achieved by considering the modified observables TO (Kuvv) + TO»(Kgyy) and TO (Kuvv)— TO»(Kgvyvy) 
in the current analysis. 


The analysis is performed in several steps. First, multi-dimensional histograms of observables are created 
in 81 bins of Kyyy/Ksm and (KAyy/Ksw) : tan a for all fits. The predicted shapes of the observables for the 
signal are produced by reweighting the base Monte Carlo sample described in Section 4. The correspond- 
ing weights are derived using the analytical calculation of the H — ZZ* — 4€ matrix elements at leading 
order in perturbative QCD. The weights are calculated and applied at the Monte Carlo generator level. 
The observables used in the analysis are evaluated after detector simulation, accounting for the detector 
acceptance, resolution and reconstruction efficiency. The distributions of observables for backgrounds 
are estimated using Monte Carlo (for the irreducible background) and data-driven techniques (for the 
reducible backgrounds) described in Section 5 and Refs. [12, 18]. 


The distributions of observables are three-dimensional: TO (Kayy, œ), TO2(Kayy, à), BDT(ZZ) and T O ((Kriyv) 
+ TO2(Kyvv), TO i(kuvv) — TO»(Kgvy), BDT(ZZ) respectively. To obtain a reliable description for bins 

with an insufficient number of Monte Carlo events, the Kernel Density Estimation [43] smoothing pro- 
cedure is applied to signal and background multi-dimensional histograms. In the smoothing procedure 

the smearing is done separately in four bins of BDT(ZZ), preserving the original normalisation. 


24 


€ A4B5ET "rr pe Tp UTQ[UTTT[TTTTTTTTTTAN 
Iz 


o. p ] m 
s 35r (s = 7 TeV, 4.5 fb ATLAS + o E js = 7 TeV, 4.5 fb" ATLAS J 
E 30. (s = 8 TeV, 20.3 fb" HoOZZ*4l| 4 $ 40- ls = 8 TeV, 20.3 fb” HoZZ'24 J 
£ ` _ i J E356 — : j 
W osf RE Ku! Kou = 1, Kaw 20,0 = 0 4 n cc Kuyy/Kgy = 1, Kay 20,0 = 0 J 
[ Te Kavv/ Ksm 771, Kaw 70, 01 = 0 1 90 1 Kuvy/ Ksm 771 , Kayy 70,0. = 0 E 
20- - 25- = 
15- : Bue E 
Ë | 156 4 
10- a 405 j 
5E J 5- E 


(a) (b) 
e S nea alaq Susu Sura pQ ec aT pr a QggrErtttt ptt peTUTETETTDTT T T rT 
d 25- (s = 7 TeV, 4.5 fb" ATLAS | ONE (s-7TeV, 4.5 tb" ATLAS 
@ | ís-8TeV,20.3 fb" H>ZZ* 541 J © 305 s = 8 TeV, 20.3 fb" H>ZZ*>41 4 
aol ] ai m : i 
ui 20- " Ll 25 [7 (Rayy/Kgy)  tan(a) = 5, K = 0 j 
L ] E (Kavv/ Ksm) tan(a) =5 , £ - 0 J 
18, ] 20- D : 
[ : 15- 4 
10- 4 E i 
" 1 10F 4 
5- Ë 5- E 
ok : o= : 
-3 -2 -i 0 1 2 3 -3 
TO (Rays) TOR, vy.) 
(c) (d) 
Lo E T T T T T T T T T T T T T T =| 
S ASE (s. 7 TV, 4.5 tb" ATLAS {+ 
~ 40F 4 4 * Data 
a VE Ys = 8 TeV, 20.3 fb H>ZZ*—4 4 
i 3 5E E Bl Background zZ* 
Mi E ~ Was m T Ray e 0,0 7 0 E [Bj Background Z4jets, tt 
BOE o figa ord Ray =O, nO J — SM: gy 1, Ky 0 y = 0, =0 
C 4 ` SM o’ “HVV — ' “AVV T * m 
25- E 
20: 1 
15: E 
10} j 
5 3 
0 zi 
0 0.25 0.5 0.75 1 


(e) 


Figure 8: Distributions of the observables used in the matrix-element-observable fit. (a): TO (Kuvv) + TO2(Kyvv), 
(b): TO (Kuvv) — TO2(Kyvv), for the Monte Carlo signal generated with (Kyvy/Ksm = 0, xXl;kayy = 0). (c) 
T Oi(Kavv, œ), (d): TO2(Kavv, œ) for the Monte Carlo signal generated with ((KAyy/Ksw) : tana = 0, +5; kgyy = 0). 
(e): BDT(ZZ) for the Monte Carlo signal generated with (kqyy/kKsm = 0, +1; kayy = 0). The expected background 
contributions are shown as filled histograms on each plot. 
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The final pdfs used in the fits are obtained by applying linear histogram interpolation between the multi- 
dimensional bins of Kyyy/kKsm and (Kayy/Ksm) : tana. The individual likelihood functions per centre-of- 
mass energy (+s) and final state (FS) are: 


c(a 


Sa AYU ai, J =| |” as = ESVY apa: J + bi) (11) 
KSM KSM KSM KSM 


i 
where P is the probability density function for the data vector Q, given the signal model s and back- 
ground model b. The index i runs over all the bins of multi-dimensional histograms of observables and 


Ó represents the vector of nuisance parameters corresponding to systematic uncertainties. Fits to data are 
performed by minimising the negative log-likelihood function with respect to the ratios of the couplings: 


Kuvv KAVV 4 = 
° tana, 0| = -21n | | | |z [e 
KSM KSM | 


ys FS 


k z £ 
EE AWV tan). (12) 
KSM KSM 


Lo 


The test statistic g^ = —21n(A) is defined as the profiled value of L of Eq. (12). To ensure the correctness 
of the statistical treatment and the absence of significant biases, a series of tests were performed before 
applying the fit to the data. Asimov datasets [31, 32] created from independently generated Monte Carlo 
samples with Kyyy/Ksm and (Kayy/Ksm) : tana equal to 0, +2, +4, +6, +8 and +10 were injected into the 
analysis procedure. The tests were repeated for samples corresponding to 1 and 100 times the LHC Run-I 
integrated luminosity. In all cases the fitted values of coupling constants were found to be in agreement 
with the injected values within statistical uncertainties. 


The results of the matrix-element-observable fit were validated and cross-checked using a nine-dimensional 
matrix-element method (9D fit). The method implements a multivariate per-event extended likelihood that 
is sensitive to both the Kyyy/ksm and (Kavyy/KksM) : tan o mixing parameters and is based on nine exper- 
imental observables. The probability model is constructed with separate components for signal, the SM 
ZZ* background and the reducible background. The background components are assumed to be inde- 
pendent of the Higgs boson tensor structure, so all of the sensitivity to mixing parameters comes from the 
signal component. Each component depends on nine experimental observables: mae, prr, Nae, Cos 6”, 
cos 0j, cos 65, O, mio and m34 (described in Section 5.4). 


The main sources of systematic uncertainty for the tensor structure measurements are the same as dis- 
cussed in Section 5 since they are based on the same four-lepton variables. Several additional sources 
of uncertainty, specific to each of the methods, are also taken into account. For the matrix-element- 
observable fit, the uncertainty related to the Kernel Density Estimation smoothing procedure applied to 
signal and background multi-dimensional histograms is considered. To estimate the influence of this un- 
certainty on the final result, a procedure similar to the one described in Section 5 is employed. The impact 
of the different sources of systematic uncertainty on the final results is evaluated by comparing the BSM 
exclusion limits obtained with a specific systematic uncertainty included or excluded in the fit, while 
excluding all other systematic uncertainties. A similar conclusion holds in the fixed hypothesis test: the 
systematic uncertainties have a very limited impact on the final result. The most important uncertainties 
are related to the estimates of the reducible backgrounds. The relative impact of these uncertainties on the 
final 95% CL exclusion limit on BSM couplings was found to be around +1%. The second most import- 
ant group of sources of systematic uncertainty is related to the theoretical uncertainties on the production 
cross section of the ZZ* background process. Their relative impact on the final result is found to be less 
than +1%. The precision of the tensor structure analysis is thus dominated by the statistical errors. 
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In this paper, only results based on the matrix-element-observable approach are reported. The 9D ap- 
proach was used as a cross-check and produced results compatible with the matrix-element approach. 


6.4 Individual and combined results 


The results of the tensor structure analyses performed in the H —^ WW* — evyy channel are reported in 
Ref. [8] and, for completeness, they are also summarised in Table 7. 


Coupling ratio Best-fit value 95% CL Exclusion Regions 

H ^ WW* 5 evpy Observed Expected Observed 

Kuvv / KSM -1.3 [-1.2,-0.7] (—co, -2.2] U[- 1, -0.85] U[0.4, co) 
(Kavv/Ksm) - tana —0.2 n.a. (—co, —6] U[5, co) 


Table 7: Fitted values of Kgyy/ksw and (Kavy/Ksm) : tan œ and 95% CL excluded regions obtained in H —^ WW* > 
evuy analysis. The expected values are estimated for the signal strength measured in data and assuming best-fit 
values for all other nuisance parameters. Only data collected at ys = 8 TeV are used. The symbol "n.a." denotes 
the absence of 9596 CL sensitivity. 


The distributions of the test statistic for fits of Kyyy/Ksm and (Kayy/KksM) : tan œ measured in the H — 
ZZ* — 4t analysis are shown in Figure 9. The expected curves are calculated assuming the SM J? = 0* 
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Figure 9: Expected and observed distributions of the test statistic for fits of (a) Kyyv/Ksm and (b) (Kavv/Ksm) : tan œ 
for the H — ZZ* — 4€ analysis. The expected curves are calculated assuming the SM JP = 0* signal and produced 
with the SM signal strength u = 1 and with the signal strengths fitted to data. The horizontal dotted black lines 
represent the levels of —2 In 4 above which the values of coupling ratios under study are excluded above 68% and 
95% CL, respectively. 


signal, both with the SM signal strength, u = 1, and with the signal strength fitted to data, ñ. The 
fitted values of Kyyy/kKsm and (Kavy/Ksm) : tana, together with the intervals where these couplings are 
excluded at above the 95% CL, are reported in Table 8. The fitted values agree with the SM predictions 
within uncertainties. 
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Coupling ratio Best-fit value 95% CL Exclusion Regions 


H > ZZ* 4€ Observed Expected Observed 
Kuyy/KSM -0.2 (—oo, —0.75] (J[6.95, co) (—oo, —0.75] U[2.45, oo) 
(KAvy/KsM) : tana —0.8 (—eo, —2.95] U[2.95,00) (—eo, -2.85] U/[0.95, co) 


Table 8: Expected and observed best-fit values of kyyy/ksm and (Kavy/Ksm) : tan æ and 95% CL excluded regions 
obtained in the H — ZZ* — 4€ analysis. The expected values are estimated for the signal strength measured in 
data and assuming best-fit values for all other nuisance parameters. The data for ys = 7 TeV and Vs = 8 TeV are 
combined. 


The measurements from the H —^ WW* — evuv and H — ZZ* 4€ channels are combined under the 
assumption that the BSM ratios of couplings Kyyy/Ksm and (Kayy/KsM) : tan o are the same for the W 
and Z vector bosons. A common test statistic is obtained by combining the profiled likelihoods of the 
individual channels. The expected distributions of the likelihoods, for the signal strength values obtained 
from the fits to the data (u = fi), are presented in Figure 10. The observed distributions of profiled 
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Figure 10: Expected distributions of the test statistic for the combination of H — WW* — evuv and H > ZZ* — 4€ 
analyses as a function of BSM coupling ratios (a) kyyy/Ksm and (b) (Kavv/Ksm) : tana. The expected values are 
estimated for the signal strengths measured in data and assuming best-fit values for all other nuisance parameters. 
The 68% and 95% CL exclusion regions are indicated as lying above the corresponding horizontal lines. The 
individual distributions for H — WW* — evuv and H — ZZ* — 4€ channels are shown. 


likelihoods for the combination of H — WW* — evuv and H — ZZ* — 4€ measurements are presented 
in Figure 11. The asymmetric shape of the expected and observed limits in the Kyyy/Ksm results is mainly 
due to the interference between the BSM and the SM contributions that gives maximum deviation from 
the SM predictions for negative relative values of the BSM couplings. Here the signal normalisations are 
treated as independent nuisance parameters of the different decay channels and the different centre-of- 
mass energies. The other nuisance parameters related to the experimental and theoretical uncertainties 
are treated as correlated when appropriate. The resulting 9596 CL exclusion regions for the combinations 
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Figure 11: Expected and observed distributions of the test statistic for H —^ WW* — evuv and H > ZZ* — 4€ 
analyses and their combinations. The distributions are shown as a function of the BSM coupling ratios kyyy/ksm 
and (Kayy/Ksm) : tana. The 68% and 95% CL exclusion regions are indicated as lying above the corresponding 
horizontal lines. (a) and (b): individual H —^ WW* — evuv, H — ZZ* — 4t and combined observed distributions. 
(c) and (d): expected and observed combined distributions. The expected distributions are presented for the SM 
signal strength u = 1 and for the signal strengths obtained from the fit to data. 
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of H —^ WW*  evuv and H — ZZ* 4t channels are listed in Table. 9. 


Coupling ratio Best-fit value 95% CL Exclusion Regions 

Combined Observed Expected Observed 
Kuvv/KsM —0.48 (—œ, —0.55] U[4.80, œ) (—eo, —0.73] L[0.63, co) 
(Kavv/Ksm) : tana —0.68 (—eo, —2.33] U[2.30, 00) (—eo, -2.18] U/[0.83, co) 


Table 9: Expected and observed best-fit values of (a) Kgyy/ksw and (b) (Kavv/Ksm) : tan à and 95% CL excluded 
regions obtained in the combination of H — ZZ* — 4t and H —^ WW* — evpyv analyses. The expected values are 
estimated for the signal strengths measured in data and assuming best-fit values for all other nuisance parameters. 
The signal strengths are treated independently per decay channel and per collision energy. 


7 Conclusion 


Studies of the spin and parity of the observed Higgs boson in the H —^ ZZ* 4t, H —^ WW*  evyuv 
and H — yy decay processes are presented. The investigations are based on 4.5 fb ^! and 20.3 fb"! 
of pp collision data collected by the ATLAS experiment at the LHC at Vs = 7 TeV and ys = 8 TeV, 
respectively. The SM Higgs boson hypothesis, corresponding to the quantum numbers J^ = 0°, is tested 
against several alternative spin and parity models. The models considered include non-SM spin-0 and 
spin-2 models with universal and non-universal couplings to quarks and gluons. The combination of 
the three decay processes allows the exclusion of all considered non-SM spin hypotheses at a more than 
99.9% CL in favour of the SM spin-0 hypothesis. 


The tensor structure of the HVV interaction in the spin-0 hypothesis is also investigated using the H — 
ZZ* — At and H — WW* — evyv decays. Only one BSM tensor coupling is investigated at a time, while 
the other one is set to zero. The observed distributions of the variables sensitive the ratios of the BSM to 
SM tensor couplings, Kyyy/Ksm and (Kayy/Ksm) : tan œ, are compatible with the SM predictions. 


Values of the BSM tensor couplings outside of the intervals —0.75 < Kyyy/Ksm < 2.45 and -2.85 < 
(Kavv/Ksm) : tana < 0.95 are excluded at the 95% CL for the H — ZZ* — 4€ process. For the H > 
WW* — evuv process the ranges —2.2 < Kyyy/Ksm < —1.0 and —0.85 < Kyyy/Ksm < 0.4 and —6.0 < 
(Kavv/Ksm) : tana < 5.0 are excluded at the 95% CL. 


The results from the H —^ WW* — evuv and H — ZZ* — 4€ decay channels are combined under the 
assumption that the Kyyy/kKsm and (Kayv/Ksm) : tan œ couplings have the same values for the HWW and 
HZZ processes. As a result of this combination, the regions outside of —0.73 < kKyyy/Ksm < 0.63 and 
—2.18 < (Kayv/Ksm) : tana < 0.83 intervals are excluded at the 95% CL. The corresponding expected 
not-excluded intervals at the 95% CL, assuming the SM Higgs boson hypothesis and the signal strength 
values measured in data, are —0.55 < Kyyy/Ksm < 4.80 and -2.33 < (Kayy/Ksm) - tana < 2.30. 
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Appendix A 


To compare the exclusion limits obtained in this analysis to other existing studies, the final results of 
this analysis are also expressed in terms of effective cross-section fractions (fo, $72) and (74, $54). The 
definitions proposed in Section 11.4.2 of Ref. [3] and Section II of Ref. [44] are used: 


|20; : 
fa c Re to mt) as 


Igo + lgo + |g4Ë oa gi 


Here the symbols gi, g2 and g4 denote the SM, BSM CP-even and BSM CP-odd tensor couplings of 
the HVV scattering amplitude, respectively. The numeric coefficients o1, o? and c4 are effective cross 
sections of the HVV interactions calculated when only the g1-, g2- or g4-related terms are present in the 
amplitude, respectively, such that g; = 1, giz; = 0. 


When, in addition to the SM term, only one CP-even or CP-odd BSM contribution is present, the con- 
version between the parameterisation used in this analysis and the (/7;, óç;) parameterisation is given by 
Eq. (13) rewritten in the following way: 


d ri (i = 2,4) (14) 
i ; L= 2, , 
"Tod i 
where r4, and r21 are chosen such that: 
cuvy [kuvv Ceiba 
à = Cri] andy ma harf ns B 
osm À ksm osm V ksm 


The numeric coefficients osm, eg yy and c'Ayy are effective cross sections of the HVV interaction calcu- 
lated when only each of the ksm-, Kyy- and Kayy-related terms is present in the Lagrangian. 


For consistency with previous measurements reported in Ref. [5], the expected and observed results of 
the current analysis of the H — WW* — evuv and H —^ ZZ* — 4€ channels and for their combin- 
ation are expressed in terms of fç; and ¢,; parameters for the H — ZZ* — 4€ decay, ( dA) and 
( 2 ot ). These parameters are denoted hereafter by (fg2, 672) and (fg4, 64). The corresponding results 
are presented in Tables 10 and 11. To obtain these results, the effective cross sections osm, Oyyy and 
Tayy of the HZZ interaction are calculated using the MapGrapHS_aMC@NLO Monte Carlo generator 
[16] at leading order. The ratios of cross sections used in the calculation are: oyyy/Osm = 0.349 and 
Tavv/Osm = 0.143, respectively. 
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Expected 95% CL limits 
H > WW* 5 evuv 


n.a. forój -0 and fy <0.15; fj; > 0.33 for dg. = r 
n.a. for dg4=0 and n.a. for bg4 = z 
H > ZZ* > 4€ 
Je < 0.94 for dg2 — 0 and Je < 0.16 for dg. = x 
Jaa < 0.56 for 4,4 =O and Jya < 0.56 for 644 = z 
Combination of H — ZZ* — At and H —^ WW* — evyuv 
Je < 0.89 for dg2 =0 and fg» < 0.096 for g2 = x 
Jaa < 0.43 for 4,4 2 0 and Jya < 0.44 for 644 = z 


Table 10: Expected limits on (fo, $2) and (74, 6,4) parameters defined in Ref. [3] obtained in the analyses of 
the H —^ WW* — evuv and H — ZZ* — 4€ channels and for their combination. The symbol "n.a." denotes the 


absence of 9596 CL sensitivity. 


Observed 95% CL limits 
H > WW* > evuv 
Je < 9.053 for 2 =0 and fo < 0.20; 0.26 < fo < 0.63 for ó; = x 


Jaa < 0.78 for ¢g4 =O and fg4 < 0.84 for 644 = z 
H ZZ > 4€ 
Je < 0.68 for øp =0 and Je < 0.16 for 6; = z 
Jya < 0.11 forój4 =0 and Jaa < 0.54 for 644 = z 
Combination of H — ZZ* > At and H — WW*  evyuv 
Je < 0.12 for ġ2 =0 and Je < 0.16 for ó = m 
Jaa < 0.090 for øg =Ü and Jya < 0.41 for 644 = z 


Table 11: Observed imits on (f42, $52) and (fa, g4) parameters defined in Ref. [3] obtained in the analyses of the 
H 5 WW* > evuv and H — ZZ* — 4t channels and for their combination. 
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